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ABSTRACT 

The spectral and temporal properties of the non-thermal emission of the nearby XRF 060218 in 0.3-150 
keV band are studied. We show that both the spectral energy distribution and the light curve properties 
suggest the same origin of the non-thermal emission detected by Swift BAT and XRT. This event has the 
longest pulse duration and spectral lag observed to date among the known GRBs. The pulse structure 
and its energy dependence are analogous to typical GRBs. By extrapolating the observed spectral lag to 
the CGRO/BATSE bands we find that the hypothesis that this event complies with the same luminosity- 
lag relation with bright GRBs cannot be ruled out at 2a significance level. These intriguing facts, along 
with its compliance with the Amati-relation, indicate that XRF 060218 shares the similar radiation 
physics as typical GRBs. 

Subject headings: gamma-rays: bursts — method: statistical 



1. INTRODUCTION 

X-ray flashes (XRFs), cosmic explosions with lower 
spectral peak energies (E p ) in vF v spectra than typical 
gamma-ray bursts (GRBs) (Heise et al. 2001; Kippen et 
al. 2003), are thought to be the low energy extension of 
typical GRBs (Lamb et al. 2005; Sakamoto et al. 2004, 
2006; Cui et al. 2005). They may be GRB jets (uniform 
or structured) viewed at off-axis directions (e.g. Zhang et 
al. 2003; Yamazaki et al. 2004; Zhang et al. 2004) or in- 
trinsically different events (Lamb et al. 2005; Soderberg et 
al. 2005). It is long speculated that long duration GRBs 
originate from a relativistic jet emerging from a collapsing 
massive star progenitor (Woosley et al. 1993; Paczyn- 
ski 1998; MacFadyen et al. 1999; Zhang et al. 2003), 
and associations of core-collapsing supernovae (SNe) with 
long GRB afterglows have been spectroscopically identified 
in a number of systems, including GRB980425/SN1998bw 
(Galama et al. 1998), GRB030329/SN2003dh (Stanek et 
al. 2003; Hjorth et al. 2003), and GRB031203/SN 20031w 
(Malesani et al. 2004). The firm spectroscopic association 
of the nearby XRF 060218 (at z=0.0331, Mirabal et al. 
2006) detected by Swift (Campana et al. 2006) with the 
Type Ic SN 2006aj was established recently (Modjaz et al. 
2006; Pian et al. 2006; Sollerman et al. 2006; Mirabal et 
al. 2006; Cobb et al. 2006; Soderberg et al. 2006). This 
suggests that the progenitors of both GRBs and XRFs are 
related to the death of massive stars and XRFs are the 
low energy extension of the more "standard" GRBs. 

Some empirical relations have been discovered from typ- 
ical GRBs. It is interesting to verify whether XRF 060218 
satisfies these relations as an approach to access its "stan- 
dardness". It has been reported that XRF 060218 com- 
plies with the isotropic energy vs. spectral peak energy 
(Ei SO — E p ) relation derived from long GRBs (Amati et al. 



2002) and other XRFs (Amati et al. 2006). The pulses in 
GRB light curves usually display a fast-rise-exponential- 
decay (FRED) shape and it has been found that the pulse 
width is related to the energy of the observational band as 
uj oc E-° A (Fenimore et al. 1995; Norris et al. 2005). A 
correlation between the isotropic luminosity and the spec- 
tral lag (Liso — t relation) of light curves was also dis- 
covered with typical GRBs as L- mo cx r -118 (Norris et al. 
2000). Both lo — E and L- lso — r relations are related to the 
structure of light curves. Since most photons of an XRF 
are in the X-ray band, previous GRB missions did not ob- 
serve the real light curves of XRFs. This makes it difficult 
to identify the temporal structure and the spectral lag of 
XRFs. XRF 060218 is long and soft, and Swift BAT and 
XRT simultaneously collected the data in the gamma-ray 
to X-ray bands. This makes it possible to measure its tem- 
poral structure and to examine whether it complies with 
the same uj — E and L\ SQ — r relations derived from GRBs. 
In this Letter we focus on this issue. Throughout the pa- 
per Hq = 71 km s _1 Mpc -1 , Q m = 0.3, and S1a = 0.7 are 
adopted. 

2. DATA 

XRF 060218 was detected with the Swift /BAT on 2006 
February 18.149 UT. It is a long burst, with a duration 
T 9Q ~ 2000 seconds in the 15-150 keV band. Swift slewed 
autonomously to the burst and the X-ray telescope (XRT) 
and UV/Optical Telescope (UVOT) began collecting data 
159 s after the burst trigger. 

XRF 060218 was an image trigger. The BAT event 
data lasted only until t ~ 300 seconds after the trigger. 
The survey data are used to derive the BAT light curve. 
XRF 060218 is very soft, with most of the emission in 
the BAT band being lower than 50 keV (Campana et al. 
2006). In order to obtain a high level of signal-to- noise 
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ratio of the BAT light curve, we use the light curve in 
whole BAT band, i.e., 15-150 keV. The first orbit of the 
XRT data is fully in the Windowed Timing (WT) mode. 
We extract the light curves and spectrum of the XRT 
data with the Xselect package. The spectrum is grouped 
with the tool grppha, and the spectral fitting is carried 
out with the Xspec package. Following Campana et al. 
(2006) we fit the XRT spectrum in the first orbit with 
a model combining a black body component with tem- 
perature kT and a cut-off power law (F oc E~ T e~ E / E °) 
component. The absorption in both Milky Way 1 and the 
GRB host galaxy (N^f st ) is incorporated. We obtain 



N% ost = 0.63t°;^ x 10 22 



kT 



0.122±°;°{« keV , 



r = 1.78±8;gf, and E c = 21.8l^ 5 7 8 keV (corresponding to 
E p ~ 5 keV), with reduced x 2 = 1.55 for 769 degrees of 
freedom. We derive the unabsorbed light curves in the 0.3- 
2 keV, 2-5 keV, and 5-10 keV bands from the XRT data. 
As reported by Campana et al. (2006), a thermal compo- 
nent exists, and is likely of different origin (e.g. the shock 
breakout emission associated with SN 2006aj; c.f. Li 2006; 
Ghiscllini et al. 2006b) from the non-thermal component. 
We therefore subtract the contribution of this component 
from the observed light curves. The kT and the radiation 
radii (Rbb) evolve with time during the first orbit. We 
read the values of kT and Rbb from Campana et al.(2006) 
and calculate the light curves of this component in 0.3-2, 
2-5, and 5-10 keV bands. Since the temperature of the 
black body component is below 0.2 keV, the light curves in 
the energy band higher than 2 keV is essentially not con- 
taminated by the thermal component. In the 0.3-2 keV 
band, the derived light curve of the thermal component 
continuously increase with time, which could be well fit- 
ted by log ^0.3-2 = (-10.69 ±0.09) + (0.66 ± 0.03) logt in 
cgs units. We thus subtract the contribution of this com- 
ponent from the 0.3-2 keV band light curve according to 
this fitting result. The derived non-thermal light curves in 
the three XRT bands as well as the one in the BAT band 
are shown in Fig. 1(a). To clearly view the pulse width 
and spectral lag dependences with energy, we also plot in 
Fig. 1(b) the normalized lightcurves. The characteristics of 
the light curves along with the average photon energy (E) 
of each energy band are reported in Table 1. 

3. RESULTS 

3.1. Spectral energy distribution 

One important question is whether the XRT non- 
thermal emission and the BAT emission are of the same 
origin. In order to clarify this we first study the spectral 
energy distribution (SED) of the burst using joint BAT- 
XRT data. Since this event is image trigger, the BAT event 
file contains only the data in the first 300 seconds since the 
BAT trigger. XRT began collecting data 159 seconds af- 
ter the BAT trigger. We thus only obtain simultaneous 
observations of the two instruments from 160 to 300 sec- 
onds with the event files. The joint-fit SED in this period 
is shown in Fig. 2, which is well fitted by the BB+CPL 
model (Campana et al. 2006) with the following parame- 
ters: N^ ost = 0.65t°;°^ x 10 22 cnT" 2 , kT = 0.131°;°* keV, 
1.56±g;gf and E c = 122±^° keV (corresponding to 
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freedom. The E p strongly evolves with time, from 54 keV 
at the beginning down to ~ 5 keV at later times. This is 
consistent with that reported by Campana et al. (2006) 
and Ghiscllini et al. (2006a). The model fitting results are 
shown in Fig. 2. One can observe that the BAT-component 
is a good extrapolation of the XRT non-thermal compo- 
nent. This result implies that the non-thermal emissions 
detected by XRT and BAT are of the same origin. 

3.2. Pulse width and Energy Dependence 

As shown in Fig.l, the light curves in different energy 
bands can be all modeled by a single FRED-pulse. Ko- 
cevski et al. (2003) developed an empirical expression to 
fit a FRED-like pulse, which reads, 



F{t) = F m { 



t + t 



trn + ta' L d + r d + r t 



t + to 
■t. 



(H-l)l 



E p — 54 keV), with reduced x = 0.97 for 327 degrees of 



(1) 

where t m is the time of the maximum flux (F m ), to is 
the offset time, and r and d are the rising and decaying 
power-law indices, respectively. We fit the light curves 
with Eq.(l) and then measure the pulse width, rising and 
decaying times at the full-width half-maximum (FWHM) 
of the fitting light curves, and the rising-to-decaying time 
ratio (if). The errors of these quantities are derived from 
simulations by assuming a normal distribution of the er- 
rors of the fitting parameters. The reported errors are at 
lcr confidence level. The results arc tabulated in Table 
1. We show w as a function of E in Fig. 3 (a). Appar- 
ently the two quantities are correlated. A best fit yields 
ui oc £;-°-3i±o.o3 rp^g p p arame ter ranges from 0.43 to 
0.59. It is found that XRF 060218 roughly satisfies the 
same u> — E relation (Fenimore et al. 1995; Norris et al. 
2005) and its <p values are also well consistent with that 
observed in typical GRBs (e.g. Norris et al. 1996; Liang et 
al. 2002), although it has a much longer pulse width than 
other single-pulse GRBs. These results imply that XRF 
060218 may be an extension of GRBs to the extremely 
long and soft regime. 

3.3. Spectral Lag and Energy Dependence 

The light curves shown in Fig. 1(a) display a significant 
spectral lag (r) , with soft photons lagging behind the hard 
photons, as usually seen in long GRBs (Norris et al. 2000; 
Yi et al. 2006). We illustrate this lag behavior with the 
intensity- normalized light curves in Fig. 1(b). The light 
curves peak at 405 ± 25, 735 ± 9, 919 ± 7, and 1082 ± 13 
seconds, respectively, in a sequence of high energy band to 
low energy band as shown in Fig. 3(b). The best fit to the 
correlation between the peak time (t pea k) and the average 
photon energy yields, 

log t peak = (3.04 ± 0.04) - (0.25 ± 0.05) log E. (2) 

A simple estimate of the lags between any pairs of the 4 
light curves obtains T pea k — 163 ~ 677 seconds, being con- 
sistent with that shown in Gehrels et al. (2007). We note 
XRF 060218 becomes the new record-holder of the long- 
lag, wide-pulse GRBs. The previous record holder was 
GRB 971208, with r - 58 seconds and u = 395 seconds 
(Norris et al. 2005). 



lr The Njj of our Galaxy (Dickey ct al. 1990) for this burst is ~ 1.1 X 10 21 cm- 2 
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We also calculate the lags with the cross correlation 
function (CCF) method. The errors of lags are evalu- 
ated by simulations. The results are also reported in 
Table 1. The lag derived by the CCF method (t CC f) 
is strongly correlated with r pca k, but is systematically 
lower 2 than T pea k [Fig. 4(a)]. A best fit gives tqcf = 
(-100±17) + (0.91±0.08)r pcak . 

The L lso — t relation was discovered with six bright 
BATSE GRBs (Norris et al. 2000) and the spectral lag 
was defined by the light curves in the 25-50 keV and 100- 
300 keV bands. We investigate whether the lag behavior 
of XRF 060218 is consistent with the L- iso — r relation. 
Since XRF 060218 is a soft XRF and the emission in the 
100-300 keV band is too weak to derive a light curve, we 
assume that t pea k of the light curve in the 100-300 keV 
band follows the t pea k — E relation (Eq.2) and perform the 
extrapolation. With the extrapolated t pea k we then esti- 
mate Tp ea k for the light curves in the 25-50 keV (average 
energy 30 keV) and 100-300 keV (average energy 200 keV) 
bands. We obtain r pea k = 177 ± 16 seconds. Since tccf is 
more reliable, we use the r poa k — tccf relation [Fig. 4(a)] 
to derive tqcf = 61 ± 26 seconds. This lag is used in the 
L — t relation analysis. Using the peak fluxes in the BAT 
and XRT band, we estimate Lj so = 1.2 x 10 47 ergs s _1 . 

Figure 4(b) shows the Li SO — r relation derived by Nor- 
ris et al. (2000) compared against XRF 060218 as well as 
two other nearby GRBs, 980425 and 031203. The data 
of the previous GRBs are taken from Norris et al. (2000) 
and Sazonov et al. (2004) 3 . The grey band and the two 
dashed lines mark the best fits at the la and 2a confidence 
level, respectively, and the solid line is the regression line 
for the 6 GRBs that were used to draw the L — r correla- 
tion, i.e. logi^ = (50.22 ± 0.32) - (1.21 ± 0.21) x logr 
(errors are at the la level). We can see that XRF 060218 
is definitely inside the 2a region and is marginally at the 
la region boundary. Therefore, the hypothesis that XRF 
060218 follows the L — t relation cannot be ruled out at 
the 2a significance level. We caution that the r is inferred 
from the extrapolation of the i poa k — E relation. This in- 
troduces uncertainties in deriving the lag. The other two 
nearby GRBs, 980425 and 031203, arc out of the 2a region, 
which are identified as significant outliers of this relation 
(e.g. Sazonov et al. 2004). 

4. CONCLUSIONS AND DISCUSSION 

We have investigated the non-thermal emission of XRF 
060218. The early SED of this event from 0.3-150 keV 
observed by BAT and XRT suggests that the non-thermal 
emission detected by the two instruments are the same 
component. By subtracting the contribution of the ther- 
mal emission we derive the light curves of the non-thermal 
emission. They are composed of a broad single pulse, and 



the energy dependences of the widths and the rising-to- 
decaying-time ratio of the pulses are roughly consistent 
with those derived in typical GRBs. The light curves 
show significant spectral lags, with a well-defined peak 
time sequence from high energy band to low energy bands, 
i-e. t P oak oc E~°- 25±0 - 05 . We infer the spectral lag in the 
BATSE bands and find that the hypothesis that this event 
complies with the L iso — r relation with typical GRBs can- 
not be ruled out at the 2a significance level. 

These intriguing facts, along with its compliance with 
the Amati-relation, strongly suggest that GRB 060218 is 
a "standard" burst at the very faint, long, and soft end 
of the GRB distribution. Since all these relations concern 
the temporal and spectral properties of emission, they are 
likely related to the radiation mechanisms. The results 
therefore imply that XRF 060218 and other XRFs may 
share the similar radiation physics (e.g. synchrotron or 
inverse Compton scattering in internal shocks, Meszaros, 
2002; Zhang & Meszaros 2004; Piran 2005; Meszaros 2006) 
with harder GRBs . 

As discovered by Norris (2002), the proportion of long- 
lag bursts within long-duration bursts increases from neg- 
ligible among bright BATSE bursts to ~ 50% at the trigger 
threshold, and their peak fluxes are ~ 2 orders of magni- 
tude lower than those of the brightest bursts. This argues 
that they are intrinsically under-luminous. Taken together 
with the fact that three nearby GRBs, 980425, 031203, and 
060218, are long-lagged and under- luminous, an intuitive 
speculation is that long-lag bursts are probably relatively 
nearby (e.g., Norris et al. 2005). The local GRB rate 
of these GRBs thus should be much higher than that ex- 
pected from the high luminosity GRBs (Liang et al. 2006; 
see also Cobb et al. 2006; Pian et al. 2006; Soderberg et 
al. 2006). A possible scenario to explain their wide-pulse, 
long-lag, and under- luminous features is the off-axis view- 
ing angle effect (e.g. Nakamura 1999; Salmonson 2000; 
Ioka & Nakamura 2001). Another scenario is that these 
features are intrinsic, being due to their lower Lorentz fac- 
tors (Kulkarni et al. 1998; Woosley & MacFadyen 1999; 
Salmonson 2000; Dai et al. 2006; Wang et al. 2006). 
They might be from a unique GRB population (Liang et 
al. 2006) having a different type of central engine (e.g. 
neutron stars rather than black holes) from bright GRBs 
(e.g. Mazzali et al. 2006; Soderberg et al. 2006). 

We thank the anonymous referee for helpful sugges- 
tions, and S. Campana, D. Burrows, J. Nousck, K. 
Page, T. Sakamoto, X.-Y. Wang, and Z. Li for dis- 
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2 We here derive tccf from the peak of the CCF without considering the side lobe contribution of the CCF. A fit to the CCF with a cube 
or quartic function gives a larger lag by considering the side lobe contribution (Norris et al. 2000). However, this method strongly depends on 
the artificially-selected range of CCF for the fitting. Since the light curves are a smooth pulse and their lags are significantly larger than the 
time bin, the peaks of CCFs are robust to estimate the lags. 

3 The peak photon fluxed in the 50-300 keV bands of these GRBs are converted to energy fluxes with the spectral index presented in Friedman 
& Bloom (2005), and their Lj so 's are recalculated with the cosmological parameters used in this paper. 
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Table 1 

Temporal structures and spectral lags of the light curves 



Temporal Structures 



Band 


Peak(s) 


W (s) 


Rising time (s) 


Decaying time (s) 


f 


E (kcV) 


(1)15-150 kcV 


405(25) 


889 (244) 


311 (28) 


578(185) 


0.54(0.18) 


36.9 


(2)5-10 keV 


735(9) 


1278 (45) 


475 (12) 


803(35) 


0.59(0.03) 


6.9 


(3)2-5 keV 


919(7) 


1707 (40) 


624 (8) 


1084(34) 


0.58(0.02) 


3.1 


(4)0.3-2 keV 


1082(13) 


2625 (125) 


794 (14) 


1831 (112) 


0.43(0.03) 


0.7 



Spectral Lags 



Bands 


A£(kcV) 


Tpcak(s) 


tccf(s) 


Bands 


A£(keV) 


Tpeak(s) 


tccf(s) 


(l)-(2) 

(1) -(4) 

(2) -(4) 


30.0 
36.2 
6.2 


330(26) 
677(28) 
347(16) 


249(37) 
518(70) 
173(25) 


(1) -(3) 

(2) -(3) 

(3) -(4) 


33.8 

3.8 

2.4 


514(26) 
184(11) 
163(15) 


389(47) 

81(12) 

43(8) 



Note: all the errors are derived by simulations and in la significance level. 
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Fig. 1. — (a) Unabsorbed light curves of the non-thermal gamma-rays/X-rays in the energy bands of 15-150 keV, 5-10 keV, 2-5 keV, and 
0.3-2 keV, respectively. The fitting curves with Eq. (1) are plotted, (b) Normalized light curves from the empirical model fitting. 
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Fig. 2.— The BAT-XRT joint spectral energy distribution from 160 to 300 seconds since the BAT trigger. The BB+CutoffPL fitting model 
is also shown. 
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Fig. 3. — The pulse duration (panel a) and the peak time (panel b) as a function of the average photon energy of the non-thermal emission. 
The solid lines in both panels are the best fits. 



9 




Fig. 4. — Panel (a): Comparison of the spectral lags derived from the peak times and from the CCF method. The solid line is the best fit. 
Panel (b): Isotropic gamma-ray luminosity as a function of spectral lag. The spectral lags of typical GRBs and GRB 980425 are calculated 
with the light curves in the 25-50 keV and 100-300 keV bands observed by CGRO/BATSE. The lag of GRB 031203 is calculated with the 
light curves in the 20-50 keV and 100-200 keV bands. The grey band and the two dashed lines mark the best fits at the la and 2a confidence 
level, respectively, and the solid line is the regression line for the six typical GRBs presented in Norris et al. (2000). 



